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We study an attractive version of the East model using the real-space renormalization group 
(RG) introduced by Stella et al. The former is a kinetically constrained model with an Ising-like 
interaction between excitations, and shows striking agreement with the phenomonology of attractive 
colloidal systems. We find that the RG predicts two nonuniversal dynamic exponents, which suggests 
that in the out-of-equilibrium regime the model coarsens via a two-stage mechanism. We explain 
this mechanism physically, and verify this prediction numerically. In addition, we predict that the 
characteristic relaxation time of the model is a non-monotonic function of attraction strength, again 
in agreement with numerical results. 
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I. INTRODUCTION 



In this paper we study the attractive East model in- 
troduced by Geissler and Reichman [lj , a kinetically con- 
strained model designed to capture the physics of col- 
loidal suspensions with short-ranged attractions 0]. Ki- 
netically constrained models (KCMs) [1 are sys- 
tems in which certain trajectories between configurations 
are suppressed 0]. As a result, they possess interest- 
ing slow dynamics H U El El El El 0. Simple 
KCMs, such as the facilitated spin models introduced by 
Fredrickson and Andersen [3j (hereafter the FA model) 
and Jackie and Eisinger |{| (hereafter the East model) 
display the slow, cooperative relaxation characteristic of 
supercooled liquids near the glass transition lR3j,|l4l . For 
a review of the glassy dynamics of KCMs see E3 • 

The attractive East model possesses in addition to a 
kinetic constraint a static attraction between facilitating 
excitations. It was shown numerically that this model 
captures some key features of the phenomenology of at- 
tractive colloids, such as a re-entrant (non-monotonic) 
timescale associated with relaxation in equilibrium, and 
logarithmic relaxation near re-entrance Q- 

In this paper we study the attractive East model using 
the real-space RG scheme of References [2I • Our rea- 
sons for doing so are twofold. First, we wish to determine 
analytically the nature of the model's out-of-equilibrium 
coarsening behaviour. We expect this behaviour to car- 
icature that of a real colloid whose attraction strength 
is fixed but whose equilibrium free volume alters sud- 
denly in response to a pressure change. Such experi- 
ments are needed to distingush between various scenar- 
ios of attractive glassiness that account for equilibrium 
behaviour. The equilibrium behaviour of the attractive 
East model, modeling a set of experiments in which free 
volume remains fixed but attraction strength varies, was 
explored in Reference pj. Second, we wish to test the 
utility of the real-space RG scheme, which has previously 
been applied to quantum spin systems [2^, reaction- 
diffusion systems |20j and simple kinetically constrained 



models E3- 

Our key results, which we present in Section I VII are 
as follows. We show that the first-order RG scheme pre- 
dicts two nonuniversal (parameter-dependent) dynamic 
exponents, which suggests that in the out-of-equilibrium 
regime the model coarsens via a two-stage mechanism. 
On short wavelengths, relaxation is approximately dif- 
fusive; on longer wavelengths, relaxation is dramatically 
sub-diffusive. We explain this mechanism in terms of a 
competition between the static attraction and the kinetic 
constraint, and verify this prediction numerically, as il- 
lustrated in Figure |21 In addition, using an uncontrolled 
second-order adaptation of the RG scheme, we show in 
Section lYlII that the characteristic relaxation time of the 
model is a non-monotonic function of attraction strength, 
in agreement with the numerical results of Ref . 

This paper is organised as follows. In Section [H] we de- 
fine the attractive East model, and for convenience cast 
its evolution operator in the guise of a quantum spin 
model. In Section ITTT1 we explain the RG scheme of Ref- 
erences [2I • In Section IIVI we apply this RG to the 
ordinary East model, emphasising the physical interpre- 
tation of the scheme, and in Section[V]we do the same to 
its attractive counterpart. We thus demonstate that the 
formalism can, with little modification, treat both mod- 
els. Those readers interested principally in the results of 
our study should focus on Sections IVII and fVIII in which 
we present our findings, to first- and second-order respec- 
tively, summarised above. We conclude in Section lVlllI 



II. THE MODEL 

The attractive East model introduced in Reference pj 
is a chain of N binary occupancies Uj = 0, 1 with periodic 
boundary conditions, an interacting reduced Hamiltonian 



N N 

H = -e riini + i + h(e, c) rn 



i=l 



(1) 



where e > 0, and the kinetic constraint of the East 
model 5]. This constraint dictates that a cell i may 
change state only if it possesses a left neighbour in the 
excited state (n.j_i = 1). We consider these occupancies 
to represent immobile (rii = 0) or mobile (rii = 1) re- 
gions of the colloid, defined over a suitable microscopic 
coarse-graining time 0, 0] . We define the equilibrium 
concentration of mobile regions as c, and we shall refer 
to occupancies = 1 as 'excitations'. 

The attraction e induces correlations between excita- 
tions. It should not be taken to represent a literal attrac- 
tion between mobile regions of a real colloid, but instead 
to describe the tendency of colloidal particles to attract 
one another ^| . These attractions may be mediated by 
dissolving a polymer in the colloid, or by changing the 
concentration of a dissolved salt solution. We regard the 
field h(e, c) as an auxiliary variable that depends on c 
and e in such a way that c is unchanged by varying e: in 
other words, varying the strength of the inter-molecular 
attraction does not change the system's free volume. The 
required adjustment of field with attraction strength fol- 
lows from the standard result for the magnetization of an 
Ising model in d = 1 i 



h(e, c) = e + 2 sinh 1 



1 - 2c 



e/2 



(2) 



The dynamics of the model is governed by a master equa- 
tion 

d t P(n, t) = -J2 w(m)P(n, t)+J2 w(l - rn)P(n', t), 

i i 

(3) 

where P(n, t) is the probability that the system has 
configuration n = {ni, . . . , m, . . . , n^} at time t (n' is 
the configuration n with occupancy rii — > 1 — n^), and 
w(rii) = w(rii, {rij}) is the probability per unit time that 
rii will change state. The {rij} are the nearest neighbours 
of i. 

We now pass to a quantum formalism |19| in the man- 
ner described in 0. Equation © is recast as the Eu- 
clidean Schrodinger equation 



(4) 



where the Liouvillian Ci is 

£, = Ci({ Sj }) | e -K«(»')-H(n)) _ e - !(«(„)-«(„<)) £.j )(5) 

and &i is the spin- flip operator defined via <Jif(n) = 
/(n')<Ti. The state vector \^(t)) can be written as 



(6) 



with |n) = The kinetic constraint is 

Ci({rij}) — rii-i, and so suppresses the dynamics of cell 
i if cell i — 1 is unexcited. 



We can write the Liouvillian J5J as 0] 

d = li ® • • • (g) (g) ^^( 2 ) ® (g) l i+2 (g> • • • <g> ljv 
+li (g) • • • (gi rij-i g 4,(1) g> g l 4+2 g) ■ • ■ <g> 1jv,(7) 

where Vi = 1 — rij is the vacancy operator. The full 
Liouvillian is composed of a sum of terms like those in 
0, with one factor for each site of the lattice. All factors 
in J7J) are 2x2 matrices. 

We can write the Liouvillian schematically as 



(8) 



where we have dropped site labels for brevity. Using the 
representation 



(9) 
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we can write the single-site Liouvillians as 

p h/2-kt/2 _ e -h/2+ke/2 



h/2-ke/2 -h/2+ke/2 I , (io) 



where the excitation number k = 1,2. The off-diagonal 
elements of Ik control the rates for the flipping processes: 
if rii has two neighbouring excitations [k = 2) it will 
excite with rate e _ ' l / 2+e , and de-excite with rate e h / 2 ~ e 
[note the global minus sign in Equation |0J]. If only its 
left neighbour is excited (fc = 1), rii will excite with rate 
e -/i/2+e/2 anc j d e _ exc jt e w ith rate e h ^ 2 ~ e ^ 2 . Since h— e is a 
monotonically decreasing function of e , the attraction 
e encourages excitations to congregate. 

The diagonal elements of £k are such that the sum of 
elements in each column is zero, a property that is re- 
quired by any probability-conserving stochastic process. 
We shall call this property stochasticity |20j , and we shall 
require that it be preserved under renormalization. 



III. REAL-SPACE RENORMALIZATION 
GROUP 

We shall study the attractive East model using a real- 
space RG developed in the 1980s, and since applied to 
quantum spin models , reaction-diffusion systems [2(j 
and nonequilibrium exclusion models |2lj |. The proce- 
dure is as follows. First, one partitions the lattice into 
blocks of b spins; we shall take 6 = 2. Then one splits 
the Liouvillian C into a reference, intra- block piece Co, 
and an inter-block interaction V in which an expansion 
will be made. By performing a suitable coarse-graining 
from the original lattice to a renormalizcd block-spin lat- 
tice, one obtains a coarse-grained or renormalizcd Liou- 
villian from which one may infer the scaling properties 
for the model in question. In Ref. 0] we showed that 
this scheme yields the low-temperature critical behaviour 
of some simple kinetically constrained models in one di- 
mension. 
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Because the attractive East model possesses 3-site in- 
teraction terms we must go to 0(V 2 ) in the perturbation 
expansion. We shall outline the procedure; see [23 for 
more details. 

We first split the Liouvillian into a reference piece and 
a perturbation, 

C = Co + V. (11) 

This splitting is arbitrary, and we shall choose the refer- 
ence piece Cq to be the simplest term that we expect to 
capture the physics of the model under study, subject to 
the requirement that it consist only of intra-block opera- 
tors. The interaction V, which constitutes the remainder 
of the Liouvillian, includes in addition inter-block inter- 
actions. For example, Cq might look like 

£ = (n®'0<8(l®l), (12) 

where the brackets indicate the way in which we choose 
to partition sites into renormalized cells. The inter-block 
interaction V might then look like 

V = (1 ® n) ® {£ ® 1) , (13) 

where now the nontrivial operators n and £ are split be- 
tween the notional renormalized cells. 

Next, we introduce the right eigenstates \<f>i) = 
ildj) i l e ')} °f t ne reference Liouvillian, 

Co \<f>i) = Efr fa) , (14) 

where \<f>i) includes both ground states \gj) with E gj = 0, 
and excited states \ei) with E ei ^ 0. We shall find 
that the reference state we choose for the attractive East 
model has three ground states (j = 1,2,3) and one ex- 
cited state (I = 1). Ground states describe slow pro- 
cesses, and excited states describe fast processes. Note 
that here and subsequently the excited states to which 
we refer are excited states of Liouvillians (evolution op- 
erators), rather than Hamiltonians (energy functions). 
The original RG scheme [2^ indeed involved Hamilto- 
nians, whereas the subsequent extension of the scheme 
to nonequilibrium processes (2jj, used in this paper, ex- 
changed Hamiltonians for Liouvillians by exploiting the 
similarity between equations such as 10} , and bone fide 
Schrodinger equations. 

We require also the left eigenstates of the reference 
Liouvillian, 

(<j> i \C = ((j> i \E <f>i , (15) 

which are generally distinct from the right eigenstates for 
non-Hermitian Liouvillians. 

The renormalization is executed by projecting the Li- 
ouvillian onto an arbitrarily-chosen coarse-graining sub- 
space S of the reference piece, corresponding to its 
ground- or low-lying excited states |20l l22l |. What re- 
mains is an evolution operator for slow, coarse-grained 
dynamical processes; whether this captures successfully 



the physics of the model depends on how valid are the 
arbitrary choices one makes, both of the reference piece 
and of the coarse-graining subspace. 

We shall construct the coarse-graining subspace S from 
only the ground states of the reference Liouvillian, in 
which case we have [2(3, [2^ 

C' (i,j) = 0; (16) 
C[(i,j) = (GilVlGj); (17) 

4M = -^(GilVltJ-^ifalVlGj), (18) 
k ^ 

where the \G a ) (and their left counterparts) are linear 
combinations of the ground-state vectors of Co . The sub- 
scripts on the matrix elements C(i,j) refer to the order of 
the perturbation series. The second order result p8|) is 
required only if the Liouvillian of the model under study 
possesses n-site interaction terms, with n > 2, as is the 
case for the attractive East model. For the ordinary East 
model, which contains only two-site interactions, the first 
order term l(T7|) is sufficient. The sum in (|T%)l runs over 
all eigenstates of Co not assigned to the subspace iS. We 
shall call this the complementary subspace, S. 

The factor of in JTSJl weights the relative contri- 
butions of excited states to the renormalization proce- 
dure. States with small eigenvalues, those corresponding 
to slower processes, are accorded more importance as a 
consequence of this factor. 

The excited states divide into two classes: those con- 
taining single-cell excitations, and those containing two- 
cell excitations. Higher-order contributions vanish by 
virtue of the orthogonality between distinct eigenvectors, 
and the fact that the model we study contains at most 
three-site, or two-cell interactions. To second order we 
acquire contributions from one-cell and two-cell excita- 
tions. 

Note that according to the RG prescription of Stella ct 
al., the sum in the second-order result (|18|) runs over all 
eigenstates of the reference Liouvillian Co not assigned to 
the subspace S. In our treatment of the attractive East 
model we shall arbitrarily discard one of the right ground 
eigenstates of Cq, I53), and we shall explain in detail in 
the following section our physical reasons for doing so. To 
first order in the RG scheme this omission introduces no 
mathematical complications. However, to second-order, 
the sum in Equation (|18l) becomes ill-defined. We find 
that Co admits one eigenstate with nonzero eigenvalue, 
and we can incorporate this into the sum without diffi- 
culty. However, by omitting \g^) from the coarse-graining 
subspace <S, and thereby assigning it to the complemen- 
tary subspace <S, we are required by the RG to sum over 
this state in (|18fl . Two issues arise. The first is that 
while I (73) is included in the subspace S, we assign its left 
counterpart (53 1 to the coarse-graining subspace S (we 
choose to do this in order to ensure probability conserva- 
tion). Thus the projection of (53 in the complementary 
subspace S is the zero vector. The second issue is that 
the eigenvalue associated with the omitted state is zero, 
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and so the denominator E g3 vanishes. As a consequence 
of these two points, the relevant term in Equation 118|) 
contains a factor 0/0, and so is undefined. Since we can- 
not treat this term in any meaningful way, we choose to 
discard it. It should therefore be borne in mind that at 
second order our treatment is entirely uncontrolled in the 
mathematical sense. Nonetheless we shall present our re- 
sults obtained in this manner, which appear to reproduce 
the numerically-observed Q non-monotonic behaviour of 
the attractive East model's charateristic relaxation time. 

It might be possible to justify our treatment to second 
order as an embedding onto a three-state coarse-graining 
subspace £3 (spanning all of the ground space) , followed 
by a projection onto the original two-dimensional space 
S |23,|2J|. This procedure would still be uncontrolled, 
however. An alternative means of regularizing the trou- 
blesome term in Equation l|18|l might be to modify the 
physics of the model in a controllable way. The phys- 
ical root of the mathematical problem discussed above 
is the (total) asymmetry of the attractive East model's 
dynamics. One could imagine replacing the attractive 
East model with a model interpolating between it and 
an attractive version of the FA model, similar to the 
interpolation model of Buhot and Garrahan |2 5ll . The 
RG properties of the latter were studied in [l8||. where 
it was shown that the ground-state vector \g 3 ) acquires 
a nonzero eigenvalue of 0(r). Here r is an asymmetry 
parameter equal to zero in the East model limit. By 
taking the attractive East model limit of the attractive 
interpolation model in a suitable way, it might be possi- 
ble to regularize the 0/0 pathology plaguing the current 
treatment. 

In the following section we describe how this formalism 
may be used to treat the ordinary East model, empha- 
sising the physical interpretation of the procedure. In 
section [V] we apply this formalism to the attractive East 
model. 



IV. RG FOR THE EAST MODEL 

In this section we shall outline the application of the 
RG procedure described above to the East model [l8| . 
We show that it is possible to partially define the scheme 
by imposing normalization and probability conservation. 
To complete the definition we arbitrarily discard part of 
the ground-state subspace of the reference Liouvillian, in 
what follows we justify this choice physically by consid- 
ering in detail the peculiar hierarchical dynamics of the 
East model. This allows us to show that the RG for the 
attractive East model, which we discuss in the following 
section, can be determined, to first order, by straight- 
forward generalization of the result for its noninteracting 
counterpart. 

The Liouvillian of the East model is 

C = n®i, (19) 



where 

^Q:r c c )' ( 2 °) 

and c is the excitation rate of a facilitated cell. Equation 
(|T§)| describes the dynamics of a cell facilitated only by 
a left-neighbouring excitation. The RG prescription 117|) 
may be written 

C'(n')=T 1 (n',n)-C(n)-T 2 (n,n'), (21) 

where T x and T 2 are respectively 2 N ' 2 x 2 N and 2 N x 2 N I 2 - 
dimensional matrices built from the left (Ti) and right 
(T2) ground-state eigenvectors of Co- The RG prescrip- 
tion is thus a mapping from a 2 Ar -dimensional Hilbert 
space of 'real' occupancies to a 2 JV / 2 -dimensional Hilbert 
space of renormalized occupancies, which we distinguish 
with primes. 

As discussed in Ref. 0|j we nn d that T\ and T 2 have 
simple interpretations in terms of the dynamics of the 
model under study. The projection matrix T\ takes a 
'real' state {11,10,01,00} and projects it onto a renor- 
malized state. For the East model we choose that states 
with at least one excitation project onto a renormalized 
excitation. Thus T\ = tf N ^ 2 where 

il = |l'){(H| + (10| + (01|} + |0')(00| (22) 

Hence blocks of spins which can, within a single time 
step, facilitate neighbouring spins are deemed also to be 
facilitating in the coarse-grained sense. Note that if we 
were to renormalize an Ising model with no kinetic con- 
straint, then we would require that flipping a spin in an 
unrenormalized configuration would result in the flipping 
of the coarse-grained spin. This is not the case for kineti- 
cally constrained models, for which there is no symmetry 
between excitations and vacancies. 

The embedding matrix T 2 reconstitutes a real state 
from a coarse-grained state, and in a sense identifies those 
original states most important for the low temperature 
dynamics of the system. As discussed in 0], we build 
T 2 so as to respect both energetic and entropic effects. 
For the East model we choose T 2 = tf N ^ 2 , where 

*2 = T ^{A|11) + |10)}<1'| + |00)<0'|. (23) 

We have defined A = c/(l — c). The structure of 
(|23|l can be partially determined by the right ground- 
state eigenvectors of Co, together with the require- 
ments of normalization and conservation of probabil- 
ity [13, but may be simply motivated as follows. To 
find t 2 , one adds to the term |00) (0'| a sum of terms 
(l + A) _1 A" 1+n2_1 \ni, n 2 ) (l'|, where n\n 2 is an unrenor- 
malized configuration, and the power of A accounts for 
the energetic weighting of these states. Thus the state 11 
is penalized by a factor of A, wheras the state 10 receives 
no penalty. 
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According to this rule we should also include in ill'MI) 
a term (1 + A) -1 1 01) (l'|, but we choose to suppress this 
for the following reason. The East model possesses an 
hierarchical dynamics jjj Il5j| . whereby two excitations 
separated by a distance d are relaxed by establishing a 
set of isolated excitations between them, at distances d/2, 
?>d/A etc. We can incorporate this behaviour into our RG 
scheme by suppressing 'frozen' configurations 01 during 
embedding, which is permitted by the structure of the 
ground state eigenvectors of the East model reference 
Liouvillian. Under repeated application of this modified 
embedding operator we then see that the most important 
state in the dynamical sense for d = 4 is 

1011, weight 1 x c 2 , (24) 

and not, for instance, 

1001, weight x c. (25) 

The latter state is favoured thermodynamically over the 
former, because it contains one fewer excitation, but the 
latter is suppressed entropically (the factor multiplying 
c 2 ). To see this, note that in order to relax the rightmost 
excitation in state Il25[l we must excite the second spin, 
followed by the third, and then finally we may relax the 
rightmost spin. To perform a similar relaxation for state 
(|24|l we may simply relax the rightmost state. Thus in 
an approximate sense we see that our ad-hoc choice of 
T2 respects the hierarchical dynamics of the East model. 
Note that for the FA model, which has symmetrical dy- 
namical rules, we choose to make no such suppression. 
Indeed, none is required. The ground-state subspace of 
the FA model's reference Liouvillian is only doubly de- 
generate, and these two states are natural choices for the 
two-dimensional coarse-graining subspace S. States i(?4*|) 
and (|25|l are then of similar importance to the dynamics. 

The renormalization prescription can thus be thought 
of as a scheme that picks out those dynamical trajecto- 
ries that most readily relax the system. As detailed in 
|18| we find for the East model the dynamic exponent 
z = — 1/ (In cm 2), to leading order in c, and a marginally 
unstable critical point at c = 0, at which the recursion 
relation for A = c/(l — c) satisfies A' = A + A 2 . From 
these results and the relation between excitation concen- 
tration and temperature, c = 1/(1 + e 1/,T ), the scaling 
properties of the East model at low temperature follow. 
For example, based on the existence of a critical point 
we can write down a scaling relation for the density of 
excitations n in the out-of-equilibrium regime, 

n(t,X)=e-' u h(t(l),\(£)). (26) 

Here d = 1 is the physical dimension, £ — In b parame- 
terizes the lattice rescaling parameter (now taken to be 
a real number), n is a dimensionless scaling function, 
and the arguments on the right hand side of (|26|l are 
the renormalized flowing time t(£) — e~ zi t and tempera- 
ture parameter A(^) = 1/ (A -1 — \n(£/£ )), respectively. 
Equation (|26|l relates a real system at long times and 



low densities (left hand side) to an effective, renormal- 
ized system at short times and high densities (right hand 
side). By imposing the matching condition e~ zl t = 1 we 
find in the regime <C In t <C 1/ (T 2 In 2) the anomalous 
coarsening behaviour 

n(t)~t- Tln2 , (27) 

which has been derived by other means and verified nu- 
merically [Tr| . 

V. RG FOR THE ATTRACTIVE EAST MODEL 

We now apply the formalism of Section 2 to the attrac- 
tive East model, whose Liouvillian is 

C = n®£ 2 ®n + n®£ l ®v. (28) 

In order to obtain a reference Liouvillian that contains 
only two-site interactions, we write the vacancy operator 
v in the second term on the right hand side of 12811 as 
1 — n, to get 

C = n®£ l +n®k®n, (29) 
where A = £2 — 1\. We choose as our reference Liouvillian 

£ = n<g>4- (30) 

The ground state embedding and projection operators 
are then similar to those for the ordinary East model. 
The projection operator T\ is identical, which may be 
motivated in a physical way by considering that the def- 
inition of a facilitating spin in the attractive East model 
is identical to that in its ordinary counterpart. The 
embedding operator T2 = tf N ^ 2 now accounts for the 
attraction-modified Boltzmann weights of two-site con- 
figurations: the attraction 'rewards' neighbouring excita- 
tions with an energy — 2e with respect to an excitation- 
vacancy pair. We enforce the same ad hoc suppression 
of the state 01 as we did in the case of the ordinary East 
model. We thus get 

t 2 = -^-^ (e*- h |11) + |1Q» (l'| + |00) (0'| , (31) 

where h = h(e, c). In the limit e^Owe recover Equation 
P3)l. with A — > e~ h . For the second order result (|T%1 we 

require also the excited projection matrix E\ — ef N ^ 2 
and embedding matrix E2 = ef N ^ 2 , which are deter- 
mined by 

e 1 = |l')«10|- e - e +' l <ll|), (32) 

and 

e 2 = (l + e ^)- 1 (-|ll) + |10))(l'|. (33) 
The prefactor in l|33|l ensures that e\ ■ e 2 = |1') 
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We now write our Liouvillian as C = Co + V, where 
£ = (n®h)®{l®l) (34) 

and 

V = A® B + C ® D, (35) 

with A = (1 ® n), -B = {£ x ® 1 + A <g> n), C = (n ® A) 
and D = (n (8 1). The brackets denote the blocking of 
the lattice into cells. 

We find under renormalization that to first order 

£[ = hAt 2 ®hBt 2 . (36) 

To second order, neglecting the pathological term dis- 
cussed at the end of Section ITTT1 we have 

£ 2 = hAe 2 ■ eiAt 2 ® t x Bt 2 ■ t 1 Bt 2 + 
t 1 Ae 2 ■ exCh ® t Y Bt 2 ■ t 1 Dt 2 + 
txAt 2 ■ hAt 2 ® hBe 2 ■ eiBt 2 + 

~tiAe 2 ■ e\At 2 ® t 1 Be 2 ■ e x Bt 2 + 

t 1 At 2 ®t 1 Be 2 - e 1 Ct 2 ®t 1 Dt 2 . (37) 

We can make sense of these results by noting that the 
combination t\At 2 — n' is a renormalized number op- 
erator, and t\Bt 2 = £' is a renormalized single-site Li- 
ouvillian. Hence the first order result H36|) looks like a 
renormalized (ordinary) East model, C! = n' <g> £'. The 
second order 2-site terms, the first four lines of Equation 
(|37|l . also have the form 'constraint ® rate', or n' ® £'. 
Note that the fourth term arises from a two-cell excita- 
tion, and so, by virtue of the factor of in Equation 

(|18|l . enters with a factor of ~. The final term in (|37|l is 
a three-site interaction that allows us to write C 2 in the 
form ©. 

In addition, to second order, we find two terms corre- 
sponding to operators not present in the original Liou- 
villian, namely 

C' 2 = t 1 At 2 ® t 1 Ae 2 ■ e x Bt 2 ® t x Bt 2 + 

t 1 At 2 ®t 1 Be 2 - e 1 At 2 ®t 1 Bt 2 . (38) 

For example, the second term in (|38|l represents a non- 
local facilitation of the form n' i _ 1 ® ■ ■ ■ ® £' i+1 . Hence to 
this order the renormalization procedure is not closed. 
We choose, arbitrarily, to ignore these terms, for two 
reasons. First, in order to treat them properly within 
the framework of the RG we would have to insert them 
into the unrenormalized Liouvillian, which, because of 
the terms' nonlocal nature we deem to be unphysical; 
and second, omitting these terms seems not to adversely 
affect our results. 



VI. RESULTS OF OUR STUDY TO FIRST 
ORDER 

With the renormalized Liouvillian now in hand we can 
extract the scaling properties of the attractive East model 



in the limit of small defect concentrations and long times. 
The first order result (|36[l shows that in the long time- 
and large length-scale limit the model behaves like an or- 
dinary East model with renormalized parameters. This 
agrees with the analysis of , which concluded that the 
attractive East model behaves on long length scales as 
an ordinary East model, albeit with a rescaled typical 
width of excited domains. We shall use the second order 
RG result l|37[) to infer the 'flow' of the model towards 
this 'renormalized' East-like behaviour. To extract the 
dynamic exponent, however, it is sufficient to use the 
first-order result, which describes the slowest dynami- 
cally relevant processes. We proceed as follows. 

The dynamic exponent z, defined via t' = b~ z t, de- 
scribes the rescaling of time as a consequence of rescal- 
ing space by a factor of b . Hence we can determine z 
by studying the ratio of renormalized to unrenormalized 
rates. We define the rates r Q , a S {1,2,3,4}, for the 
four processes of the model as follows: 



Q 


Process 




1 


111 - 


■* 101 


e h/2-e 


2 


110 - 


-> 100 


e h/2-e/2 


3 


101 - 


-> 111 


e -h/2+e 


4 


100 - 


■+ 110 


e -h/2+e/2 



Let us define the renormalized rates T a extracted from 
the first order result (|36|l . Then upon rearranging — 

b~ z T we have z a — — lnr a /ln2, where r a = Ta /r a . 
The renormalized rates are cumbersome, and so we shall 
not display them explicitly. 

We find that, unlike for the ordinary East model, there 
is no common rescaling factor z. Instead we find four 
distinct z a , whose dependence upon e is dramatically 
different. Two values, z\ and Z4, we discard as unphys- 
ical, being negative and tending asymptotially to unity, 
respectively. The other two we plot for excitation con- 
centration c = 10 -3 in Figure ^ The larger value, Z3, is 
non-monotonic, and follows approximately the behaviour 
of the East-like exponent. The smaller value, z 2l tends 
to a diffusive value 2 for large attraction strengths. 

We can understand the nature of these two dynamic 
exponents from the following simple argument. The dy- 
namic exponent characterizes the typical rate t _1 at 
which a domain wall can move a given distance £, via 
t ~ l z . Let us define a diffusion constant D, such that 

f = ((r(r) - r(0)) 2 ) = 2Dt, (39) 

and r(r) is the position of a given domain wall at time r. 
Replacing the exponent 2 in Equation (|39|l by the more 
general value 2, we can write 

t~D- x £ z . (40) 

But if we consider domain wall drift as an 'activated' 
process, one inhibited by an effective barrier of size A(£), 
then we can also write 

r^r - 1 exp(AW), (41) 



G 




FIG. 1: (Colour online) Dynamic exponent z versus attrac- 
tion strength e for the attractive East model, at excitation 
concentration c = 10 -3 . The lower horizontal line marks the 
diffusive limit z = 2; the upper horizontal line is the pre- 
diction for the East model in the absence of the attraction. 
The bold dashed lines are the phenomenological predictions 22 
(lower line, green) and 23 (upper line, red). These we expect 
to be characteristic values for the dynamic exponent of the 
attractive East model when the domain sizes are respectively 
much less than and much greater than the attraction-induced 
lengthscale a(e, c). The bold blue lines are the prediction for 
2 from the first-order RG calculation. Both correctly repro- 
duduce the e — > limit, and agree with our numerical results 
(see Figure |2J . 



where ro ~ D is the rate for attempting a barrier cross- 
ing. Equating lj!U|l and l|4T]l and taking logarithms gives 



z ~ 



(42) 



Note that the effective barrier A(£) accounts for the dy- 
namics of the process under consideration. For example, 
for the ordinary East model the barrier A(^) grows log- 
arithmically with i 0- We then obtain z ~ l/(Tln2) to 
leading order in T. 

The two dynamic exponents we find in the case of the 
attractive East model arise because the attraction im- 
poses a length scale below which the East-like hierarchi- 
cal dynamics is suppressed. This follows from the fact 
that there exists in the attractive East model a charac- 
teristic width a(c,e) of excited (black) domains [lj. To 
see this, note that the energy of a configuration 10 • • • 01 
of length a is 2ft. The energy of a similar length of chain 
with all the intermediate cells excited is ft + a(h — e). 
These energies are equal when a = ft/ (ft — e). Hence a, 
which increases as e increases, sets a length scale below 
which a 10 domain wall can readily move eastwards via 
a mechanism that excites contiguous cells. See Figure El 
for an illustration of this effect. This domain wall cannot 
move freely, however; the penalty for exciting successive 
cells to the east is, from (1), ft — e. Thus domain wall 
motion on lengthscales £ < a should proceed principally 
by diffusion in a potential (ft — e)£. The largest barrier 



X 100 




f [ G . 2: An illustration of the finite thickness of excited 
(black) domain walls induced by the attraction e. We show 
here nonequilibrium space-time trajectories for the attractive 
East model, starting from excitation concentration n(0) = 0.5 
for equilibrium excitation concentration c = 0.1. We show the 
cases e — (upper panel) and e = 5 (lower panel) . Space runs 
along the vertical axis, time along the horizontal. Notice the 
emergence of a finite thickness of excited domains in the lower 
plot, for which the attraction-induced lengthscale a « 24 lat- 
tice sites. On lengthscales I < a we expect essentially diffusive 
dynamics. 



to be surmounted is therefore (ft — e)a, and so 

(ft — e)a 



Z2 



In a 



(43) 



We call this exponent Z2 as per the RG notation, and we 
expect it to dictate the dynamics of the system when the 
characteristic domain size L(t) satisfies L(t) < a. 

However, on lengthscales L(t) > a we expect relaxation 
to proceed via East-like hierarchical dynamics, suitably 
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FIG. 3: (Colour online) Continuous time Monte Carlo sim- 
ulation data for the attractive East model with c = 10 -3 , 
showing elapsed time r (vertical axis) versus mean white do- 
main length L (horizontal axis). The system size is N = 10 5 . 
We show the average of 10 runs for each attraction strength 
e — 0, 2, 6, 8, 10, 12, from left to right. The two-stage coars- 
ening mechanism is seen clearly in the three curves at largest 
e, and the crossover occurs when the mean domain length 
exceeds approximately a(c, e). Note that the domain length 
distribution broadens as the model slows, and so the mean do- 
main length becomes progressively less good as an estimate 
of the crossover. We show the value of a as vertical dotted 
lines for e = 8, 10 and 12. We show also the RG predic- 
tions for short-wavelength coarsening at large e, z — 2 (heavy 
black dashed line), and for long- wavelength coarsening for 
e = 6{z » 6.72), e = 8 (z « 7.24) and e = 10 (z w 8.37) 
(heavy red dashed lines). We show for comparison next to 
the e — 6 line the ordinary East model low-temperature pre- 
diction, 2 w 9.97 (thin blue line), which clearly disagrees with 
the data. 



rescaled to account for the attraction-induced length a. 
A characteristic value for the relevant dynamic exponent 
can be derived by considering the typical size of white 
domains in thermal equilibrium, (L(t — > oo)) = d(c, e). 
This length may be found found from the partition func- 
tion of the d = 1 Ising model. Denoting by fi the largest 
eigenvalue of the transfer matrix 




a h/2 



(44) 



we have that white domains of size L occur with prob- 
ability P(L) = cxp(— iln/i). Hence the typical white 
domain size in equilibrium is d(c,e) — 1/ln/i. 

The logarithmic barrier confronting domain walls mov- 
ing a distance d a is then A(d) — h \n(d/a)/ In 2 0, 



h\n(d/a) 
In din 2 



(45) 



We show in Figure ^ the dynamic exponents calculated 
from both phenomenological and RG predictions. In Fig- 
ure |3 we show numerical results for the attractive East 
model in the nonequilibrium regime. The coarsening 
mechanism is as we predict from RG and physical con- 
siderations: approximately diffusive relaxation on short 
wavelengths, crossing over to dramatically sub-diffusive 
behaviour on larger wavelengths. Note that the value 
of z we display for the ordinary East model is its low- 
temperature approximation. The crossover length (the 
position of the 'kink' in the plots) is consistent with 
i(t) ~ a. For large values of e the distribution of white 
bubble lengths is very broad, and the crossover begins 
while the mean domain length is still appreciably less 
than a. 

We verified also that diffusive coarsening persists 
throughout the nonequilibrium regime for systems such 
that a > d. 



VII. RESULTS OF OUR STUDY TO SECOND 
ORDER 

While the first-order result is sufficient to determine 
the dynamic exponent of the model, to first order the 
renormalized Liouvillian (|36|l looks like that of an ordi- 
nary East model, in that the rate for creating an exci- 
tation is insensitive to the state of the right neighbour. 
However, the renormalized Liouvillian to second order, 
l|H7jl. does account for the state of the right neighbour. 
We can infer from this result a characteristic timescale 
for relaxation, as well as a timescale for crossover from 
diffusive to sub-diffusive relaxation. It should be borne 
in mind that to second order our perturbation scheme is 
formally ill-defined. 

We can infer the values of the renormalized parameters 
via the matrix elements of the renormalized Liouvillian 
to second order, £' = £[ + £' 2 . In particular, we define 
the quantities 



v( 2 ) = 



r (2) 

— • 

r (2) ' IJ 

1 1 



(2) = 



_ 1 1 1 



(2) p (2) 



r (2)p(2) ■ 



(46) 



a 2e-h{c,e) 



whose unrenormalized counterparts are a = e" 
and = e _e . From the latter two expressions we solve 
for c and e as functions of a and 0, obtaining c = /(a, 0) 
and e = — ln/3. The function / is unwieldy, by virtue 
of the complicated dependence of h(c, e) upon its ar- 
guments, Equation and so we shall not display it 
explicitly. Then we define renormalized parameters via 
d = d 2 ) = /(a( 2 \/?( 2 )) and d = e( 2 ) = -ln/?( 2 ). 

In the absence of attraction the recursion relation for c 
under renormalization reads c = c/(l — c + c 2 ), encoding 
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FIG. 4: (Colour online) Negative of the loga- 

rithm of the characteristic timescale — log 10 r = 
v\\ log 10 c from the second-order RG result, for 
c = 0.01,0.02,0.04,0.06,0.08,0.10,0.12,0.14,0.16,0.18 
(bottom to top). The degree of non-monotonicity at each 
value of c is in good agreement with the numerical results 
shown in Figure 3 of Reference £jj . The difference in absolute 
value at each concentration occurs because the RG result 
refers to a generic timescale r. For any given process r g there 
exists a basic rate F(c, e) such that r 9 = T(c, e) c~"W , and 
hence the vertical offset of — log 10 r 9 relative to the generic 
value — log 10 T g is — log T(c, e). 
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FIG. 5: (Colour online) Crossover length L xo = (eg/e) 1 ^' 
calculated from the second-order RG result, for c = 10 -3 , as a 
function of e (solid blue line). We show also the phenomeno- 
logical prediction a(c = 10~ 3 ,e) (red dashed line). The RG 
result gives correctly the trend of increasing crossover length 
with increasing e, as we expect on physical grounds. However, 
additional physical input is needed to fix the absolute scale 
of the crossover. Here we obtain close agreement between the 
two curves by choosing arbitrarily e s — 10 -3 , but the RG 
result varies strongly with e s and can be made to look very 
different from a(c, e). 



an unstable critical fixed point c* = and a stable full- 
lattice fixed point c* = 1. We shall focus on the the 
regime c <C 1, and shall avoid probing the large-e regime 
where possible static critical effects intrude. 

To extract a characteristic timescale we define the 
exponent y c via d — b Vc c. From standard RG argu- 
ments [29j we then have that v± = l/y c , where v± 
controls the divergence of the correlation length £ near 
the critical point c = 0, via £ ~ c~" x . Scaling argu- 
ments [27| dictate that the characteristic timescale di- 
verges near c = as r ~ c — "H, where v\\ = zvj_. We 
plot the logarithm of this timescale, log 10 r = u\\ log 10 c, 
in Figure 0] We find a similar degree of non-monotonic 
behaviour (roughly 1 decade at the lowest concentrations 
c) to that shown in Figure 3 of Reference pj . For consis- 
tency we use the exponent Z3 calculated to second order. 
We expect this exponent to be relevant on large wave- 
lengths, for example in equilibrium for d 3> a, as was 
the case for those simulations shown in Figure 3 of Ref- 
erence 0. 

Lastly, we can infer the lengthscale at which the attrac- 
tive East model crosses over from near-diffusive to sub- 
diffusive behavour. To extract such a crossover length, 
recall that the first order RG result, Equation l|17|). looks 
like the Liouvillian of an ordinary East model (with 
rescaled parameters) , in that the rates for excitation and 
relaxation of a cell are unaffected by the state of the 
neighbour to the east. By contrast, the original Liou- 
villian ijHJ) and the second order RG result (|18J) describe 
models whose dynamics depends on the state of a cell's 
east neighbour. 

We thus identify the parameter f3 in Equation 1461) 



as a measure of the extent to which the model looks 
'East-like' (rates insensitive to the state of the east neigh- 
bour) or 'attractive East-like' (rates dependent upon the 
state of the east neighbour). We equate the former situ- 
ation with the regime of hierarchical (dramatically sub- 
diffusive) coarsening. 

In the language of RG, then, we expect sub-diffusive 
behaviour when the renormalized parameter e' becomes 
small. Let us define the exponent y e via e' ~ b~ y 'e. From 
the linearized second-order result we find e' = e/2, giv- 
ing y e = 1. It is more meaningful, however, to retain 
all orders of e. If we do this, and then iterate the RG 
until e' ~ e s , where e s is some sufficiently small value of 
the renormalized coupling e', we find the corresponding 
value of the rescaling parameter: b ~ {tg/e) 1 ^'. Since 
this occurs as a consequence of rescaling space by a fac- 
tor of b^ 1 , we infer the crossover length L xo ~ (eg/e) 1 /^. 
On lengthscales I S> L xo we thus expect hierarchical, 
sub-diffusive relaxation. We plot this crossover length in 
Figure 03 and find that for e s = we obtain good 

agreement with the phenomenological crossover length 
a(c, e). Note that since e s is arbitrary we require ad- 
ditional physical input to fix the absolute scale of the 
crossover length. The predicted trend of increasing L xo 
with increasing e, however, is illuminating. 

As an aside, and a further demonstration of the pre- 
dictive power of the RG scheme, it is interesting to note 
that a first-order real-space RG calculation for an attrac- 
tive version of the FA model suggests a 'bifurcation' at 
large e of the dynamic exponent (Figure top). Our nu- 
merics support this prediction (Figure bottom). Our 
RG treatment of the attractive FA model is a straightfor- 
ward adaptation of the procedures used in Reference 0] 
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FIG. 6: (Colour online) Top: Real-space RG prediction for 
the dynamic exponent z for an attractive FA model as a 
function of attraction strength e; c = 10 -3 . Note the bi- 
furcation for e > 8. Bottom: A test of the RG calculation 
- Elapsed time rescaled by the square of the mean domain 
length tL~ 2 versus mean domain length L for c = 10 -3 and 
e = 0, 2, 4, 6, 8, 10, 12 (top to bottom). We show averages over 
5 runs for each attraction strength, starting from initial ex- 
citation concentration n(0) = |. For the topmost curves the 
scaling prediction r ~ L 2 , for which tL~ 2 ~ const., is rea- 
sonably demonstrated. For e > 8 (the two lowest curves) we 
see a regime in which r ~ L z with z < 2, as predicted by the 
RG calculation. We show as a heavy black dashed line the 
prediction r ~ L 1 ' 60 for e = 12. 

and in this paper: we modify the embedding operator to 
account for the attraction e, and implement the RG to 
first order. 



strength. In doing so, we have shown that the real- 
space RG scheme of References |2(| ^ , subject to our 
arbitrary selection of the coarse-graining subspace, and 
an uncontrolled approximation to second order, captures 
several interesting characteristics of this model. It also 
bears out the intuitive idea of Reference [lj concerning 
the 'renormalized East-like' nature of the model at large 
wavelengths. Because the RG scheme may be motivated 
in simple physical terms, it is therefore a useful starting 
point in directing more detailed analyses, either numeri- 
cal or theoretical. 

The two-stage coarsening mechanism of the attractive 
East model is an example of two-stage relaxation induced 
by a competition between a static attraction and a kinetic 
constraint. It would be interesting to search for such a 
mechanism in real systems. The attractive East model 
was intended originally to be studied in its equilibrium 
regime, at fixed c, corresponding to fixed colloid pack- 
ing fraction. To test the results of this paper against the 
behaviour of real attractive colloids we propose the fol- 
lowing experiment: a sudden pressure increase at fixed T, 
inducing an increase in packing fraction, and hence a re- 
duction in c as the liquid equilibrates. On the basis of our 
results we would expect particle mobility to differ quali- 
tatively depending on the wavelength one probes and the 
isothermal compression line one explores 0- For weak 
attractions and for moderate attractions at large wave- 
lengths we would expect particle motion to be controlled 
principally by free volume, and thus to be sub-diffusive 
and repulsive glass-like. For moderate attractions at 
small wavelengths we would expect particles to move in 
an approximately diffusive manner, because attractions 
render the short- wavelength structure of the collid labile. 

In general terms, it would be interesting to determine 
whether two-stage relaxation like that descibed in this 
paper could be brought about by a competition between 
two static attractions, one short-ranged and repulsive, 
and the other long-ranged and attractive [30| . 

It would also be interesting to see if the two-time scal- 
ing behaviour of the attractive East model could be de- 
termined by RG or other means; the different behaviour 
exhibited by the original East model and its attractive 
counterpart in the nonequilibrium regime (Figure 3) sug- 
gests that the two-time scaling behaviour of the models 
differ. Since this regime corresponds, according to our 
mapping, to the aging regime of an attractive colloid, 
such a study would be a valuable way of comparing the 
behaviour of the attractive East model and real colloids. 



VIII. CONCLUSIONS 
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